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ABSTRACT 

We calculate the energy that baryons must inject in cold dark matter (CDM) haloes in order to remove 
centrally-divergent DM cusps on scales relevant to observations of dwarf spheroidal galaxies (dSphs). We 
estimate that the CDM haloes often associated with the Milky Way's dSphs (Mvii/M© ~ 10^"'°) require A/s ^ 
2Q53-55gj.g order to form cores on scales comparable to the luminous size of these galaxies. While supernova 
type II (SNell) explosions can in principle generate this energy, the actual contribution is limited by the low 
star formation efficiency implied by the abundance of luminous satellites. Considering that CDM's well-known 
'core/cusp' and 'missing satellite' problems place opposing demands on star formation efficiencies, existing 
observational evidences for large cores in the most luminous dSphs require that CDM models invoke some 
combination of the following: (i) efficient (of order unity) coupling of SNell energy into dark matter particles, 
(ii) star formation histories peaking at unexpectedly high redshifts (z > 6), (iii) a top-heavy stellar IMF, and/or 
(iv) substantial satellite disruption or other stochastic effects to ease the substructure abundance constraints. 
Our models show that the tension between CDM problems on small scales would increase if cored DM profiles 
were to be found in fainter dwarves. 



1. INTRODUCTION 

Cosmological N-body simulations show that if gravita- 
tional interactions between 'standard' (i.e., massive, weakly- 
interacting) cold dark matter (CDM) particles dominate struc- 
ture formation, then galaxies must be embedded in dark mat- 
ter haloes that (i) collectively follow a mass function that di- 
verges at low masses as d/V /dMvii- ^ A/yir^ (^-S- Springel 
et al. 2008 and references therein); and (ii) individu- 
ally follow mass-density profiles characterized by centrally- 
divergent 'cusps', with p{r) oc r"' at small radii (Dubinski 
& Carlberg 1991; Navarro, Frenk & White 1996, hereafter 
NEW). The dark matter haloes inferred from observations of 
real galaxies differ significantly: the estimated mass profiles 
are consistent with homogeneous-density 'cores' (Kuzio de 
Naray et al. 2008; Battaglia et al. 2008; de Blok 2010; Walker 
& Peiiarrubia 20 11; Amorisco & Evans 2012) and the satellite 
mass function is remarkably flat (Klypin et al. 1999; Moore 
et al. 1999). 

Either the dark matter is not standard CDM or, if it is, 
then non-gravitational forces must play a significant role in 
structure formation. Indeed various baryon-physical mecha- 
nisms have been proposed and demonstrated to be capable of 
fixing both problems separately: e.g. supernova explosions 
(Navarro et al. 1996, Read & Gilmore 2005; Mashchenko et 
al. 2008; Governato et al. 2010; Pontzen & Governato 2012), 
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or the orbital decay of compact baryonic objects (El-Zant et 
al. 2001; Goerdt et al. 2010; Cole et al. 201 1) can underplau- 
sible conditions flatten the central cusps of CDM-like haloes. 
Likewise, supernova feedback, inefficient cooling of interstel- 
lar media, cosmic reionization, UV radiation, and acoustic os- 
cillations have all been invoked to explain the suppression of 
galaxy formation in low-mass satellite haloes (e.g. Tassis et 
al. 2008; Bovill & Riccotti 2009; Sawala et al. 2010; Bovy & 
Dvorkin 2012). 

Notice that solutions to CDM problems on small scales 
rest upon the efficiency with which stars form in DM haloes. 
However, while core formation generally requires an efficient 
conversion of gas into stars, suppression of galaxy formation 
obviously requires the opposite. 

In light of this apparent tension, here we examine whether 
baryon-physical solutions to CDM's 'core/cusp' and 'missing 
satellites' problems are mutually compatible. For maximal 
leverage we consider explicitly the impact of the most ener- 
getically efficient baryon-driven events-supernova explosions 
of type II (SNell)'' -on the least luminous objects associated 
empirically with dark matter haloes: the Milky Way (MW)'s 
dwarf spheroidal satellites. 

2. CUSP REMOVAL AT WORK 
2.1. Models 

We proceed to derive a rough estimate of the amount 
of work involved in transforming a centrally-divergent cusp 
into a constant-density core. We adopt the following 
cosmologically-motivated halo density profile 
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° One may show that the orbital energy of compact stellar objects is orders 
of magnitude lower than associated to SNell events. 
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where po is a characteristic halo density and is a scale ra- 
dius. We assume that the original state of the DM follows a 
NFW moder, i.e. for a core radius rc = 0, and that the main ef- 
fect of baryonic feedback is to flatten the inner cusp, so that at 
r <^rc dark matter is homogeneously distributed. Clearly, this 
solution is not unique, and in principle more complex models 
(that is, with a larger number of free parameters) can be in- 
voked in order to describe the effects of baryonic feedback in 
more elaborate a detail. In practice, however, eq. (1) provides 
a reasonable fit to the cored DM density profiles found in the 
hydro-dynamical simulations of Governatoet al. (2012, G12). 
Also, one may show that the below estimates barely change 
for more complex profiles, e.g p/po ~ {r" + rj!)"'/" at r ^ r,. 

After the DM cusp is removed we assume that the halo 
settles in a new equilibrium state. Hence, we can use the 
virial theorem to derive a lower (i.e conservative) limit of 
the amount of energy required by the transformation as /S.E = 
AW /2 = (Wcore - W'cusp) /2, where W is the halo potential en- 
ergy 



W = -4ttG 



p{r)M(r)rdr, 



(2) 



and M{r) is the halo mass profile, which can be integrated 
analytically from eq. (1) 



M(r) 



X- ln(l+f/.t)+(l-2j:)ln(l+f) 

r(2+3r) 
' 2(l+f)2 



(l+m-x) 



ln(l + r)- 



(3) 



where Mo = Airpor^^, r = r/rs and x = rc/r,. 

Given that feedback mainly shapes the inner-most regions 
of DM haloes, the parameter po can be derived from eq. (3) 
under the assumption that the halo virial mass remains con- 
stant^ i.e. Mvii-xoie = Mvkxusp = Mo[ln(l +c)-c/(l +c)], where 
c = ryir/rc is the halo concentration. 

It is useful to write 



AW= ^e(;c,c). 



(4) 



Although this equation has to be solved numerically, one 
may show that Q ^ x for x ^ 1, i.e. the energy required to 
form small cores is proportional to the core size; whereas 
for X ^ I, this function approaches asymptotically the limit 
lim^^oo ~ 0.05. This indicates that unbinding the very cen- 
tral cusp is the energetically demanding part of the transfor- 
mation, whilst re-distributing DM beyond the halo scale ra- 
dius requires, in comparison, a relatively small amount of en- 
ergy, suggesting that our results will not be strongly sensitive 
to the details of the outer halo profile. 

In our models the initial state of the halo is fully specified 
by Mvir,/"vii- and c = r^h/rs- These parameters can be fixed 
by using the results from N-body cosmological simulations. 
For example, the halo virial radius rvii- = rvii-(Mvii-,z) is calcu- 
lated from Bullock et al. (2001), whereas the concentration 

^ Adopting Einasto profiles instead of NFW models has a negligible impact 
on our estimates. 

^ However, in general we expect haloes to grow in mass while forming 
stars, i.e Myii coic IZ, ^vii-.cusp, hence our estimates of A£ must be regarded as 
conservative. 
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Fig. 1. — Upper panel: Density profiles considered in this work. Thick 
solid lines correspond to an NFW model with M^-„ = 3 X IO'Mq , scale radius 
t's ~ 1.9 kpc a concentration c = 16. Red, green and blue show cored density 
profiles of haloes with the same virial mass and a core radius = 0.2, 1 .0 and 
5 kpc, respectively. Lower panel: Slope of the mass profile F = dlnM/dlnr 
as a function of radius for the above models. NFW and cored haloes have 
r < 2 and F < 3, respectively. Recent measurements of F in the Fornax 
(closed symbol, ri, ~ 713 pc) and Sculptor (open symbol, t\ ~ 226 pc) dSphs 
are indicative relatively large (rc > 1 kpc) DM cores in both galaxies (see 
WPll for details). 

c = c(Mvir,z) follows the relation proposed by Maccio et al. 
(2007)"^. Hence the amount of work required to remove the 
DM cusp (AW) is solely determined by the size of the DM 
core with respect to the original scale radius, x = Vc/rs- 

2.2. Core size 

Recently, Walker & Penarrubia (2011; WPll) have mea- 
sured the slope of the halo mass profile (F) in the two bright- 
est MW dSphs, Fornax and Sculptor (M^/Mq ~ 3 x 10^ and 
8 X 10^ respectively; de Boer et al. 2012a,b). The values, 
r = AlogM/Alogr = 2.61!|]4^ and F = 2.95!j]'^^, rule out 
NFW profiles (F < 2) at confidence levels > 96% and > 99%, 
respectively. 

We can use eq. (1) to relate the slope of the mass profile 
measured at the half-light radius (rj,) with the core size. For 
stellar components that are deeply embedded within the DM 
halo (rh ^ Ts) we have 



(5) 



Two points can be gleaned from this Equation. First, mea- 
surements of T{rh) cannot be used to derive upper limits for 
the halo core radius, as the slope approaches an asymptotic 
value F ~ 3-3/2(rh/rj) for rc ^ rj. And second, steep mass 
profiles (F > 2.5) imply that the dark matter core extends well 
beyond the luminous radius of the dwarf, i.e > r^. 

' The scatter in this relation ±0.15 dex is neglected for reasons that will 
become obvious below. 



3 



Fig. 1 illustrates the transformation defined by eq. (1). We 
adopt a fiducial NFW model for dSphs with Mvii = 3 x IO^Mq, 
rj ~ 1.9 kpc and c = 19.6 (solid line) (see Peiiarrubia et al. 
2008). Long-dashed, dotted and short-dashed lines show 
cored profiles with = 0.2, 1.0 and 5 kpc, respectively. The 
lower panel shows the slope of the mass profiles associated 
to these models. As expected from eq. (5), the values of F 
measured in Fornax and Sculptor suggest > 1 kpc. 

It is useful to define an "observationally-relevant" range of 
core sizes. Motivated by the above measurements, we set a 
minimum core size comparable to the half-light radius of the 
dwarf r^ min ~ r^. Given that stars are deeply embedded 
within the CDM halos of dSphs, ^ 0.1 (Peiiarrubia et al. 
2008), we consider minimum core sizes 70 < /"c.min/pc < 370 
for haloes with Mvii/M© ~ 10*^"". The maximum size, al- 
though poorly constrained observationally, can be safely as- 
sumed to be r^^max < /"vii - Hence, here we inspect transforma- 
tions in the range 0.1 < r^/rj < c. This uncertainty, together 
with the details of the outer halo profile (see §2.1), has a small 
relevance to our conclusions. 

2.3. Energetics of the core/cusp transformation 

Thick red lines in Fig. 2 show the range of energies re- 
quired to remove DM cusps as a function of the halo virial 
mass. As expected from eq. (4) the amount of energy re- 
quired by the transformation scales as AW oc M^;,. At a fixed 
mass the energy range (shaded area) is relatively narrow, the 
lower limit of AW being determined by the minimum core 
size, niin ^ O.lr, (see §2.2), and the upper limit being given 
by the asymptotic limit of Q for x';^ \. For ease of reference 
we also mark the case rc = l kpc with black dots. 

The CDM haloes of bright dSphs have estimated virial 
masses in the range Mvii/M© ^ 10''"'° (Peiiarrubia et al. 
2008). Thus, baryons must produce of the order of 10^^ < 
AW/(erg) < 10^^ in order to effect the profile transformation 
on the dSph mass scales. 

3. CDM PROBLEMS AND BARYONIC FEEDBACK 

3.1. Supernova energy 

The only plausible source of such an immense amount of 
energy appears to be the averaged Esh ~ 10^' erg of kinetic 
energy released per SNII explosion (Utrobin &. Chugai 201 1). 
The fraction that contributes to the pressure term in the hydro- 
dynamical equations is the so-called 'energy coupling' (esn)- 
Its value is still a matter or debate. Typical values range 
from esN = 0.01 (e.g. Kellermann 1989) to 0.40 (Governato 
et al. 2010). Recently, Revaz & Jablonka (2012) argue that 
esN ~ 0.05 in order to describe the relation between metallic- 
ity and luminosity found in dSphs (e.g. Kirby et al. 2011). 
Note, however, that the fraction of energy transferred into the 
gas (esn) may be much larger than the fraction of energy even- 
tually transferred into the DM distribution, which we term 
cdm- We will adopt epivi = 0.4, which gives an upper limit to 
the energy available, since edm < esn, and a value esN = 0.4 is 
already quite optimistic (but facilitates comparison with Gov- 
ernato et al results). 

Notice also that infening the presence of tiny (Tc /"h) DM cores in 
dSphs is extremely challenging. 
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Fig. 2. — Minimum energy (Ail, red dashed ai'ea) required to form a DM 
core with a size 0.1 < rr/r, < c, where c = r^\Jrs is the halo concentration 
an NFW halo with virial mass Myir at z, = 0. The energy required to form 
a DM core with = 1 kpc is marked with black dots for reference. On 
the right vertical axis we plot the luminosity derived from a star formation 
efficiency iv = iv(Mvi, ) tuned to reproduce the number of luminous satellites 
in our Galaxy. Luminosities are converted into SNell energy using eq. (6) 
and adopting a strong energy coupling eoM = 0.4. Notice that the tension 
between the 'core/cusp' and 'missing satelHte' problems arise in haloes with 
Mvir < lO'^M©. 

The stellar mass (M*) required to generate an energy /S.E 
through SNell explosions is 

A£ = — > 8Mo)£sNeDM; (6) 

where we assume that stars form following a universal Initial 
Mass Function (IMF in short), ^(m*); and that only stars with 
masses m^, > 8Mq undergo SNell during their last evolution- 
ary stages. For simplicity we adopt a Kroupa (2002) IMF, 
which gives a fraction of massive single stars ^(m^ > 8M0) = 
0.0037 and a mean stellar mass (m,t) = 0.4Mq. 

Eq. (6) relates stellar mass to the amount of SNell energy 
available to do work on the halo, and we use it to draw the 
right-vertical axis in Fig. 2. If we allow stellar mass fractions 
to take values as large as the universal baryon fraction -i.e. 
M*/Mvir = nh/n,„ ~ 0.047/0.28 ~ 0. 16-then Fig. 2 indicates 
that core formation via SNell feedback is indeed energetically 
viable on all mass scales. However, this quantity is strongly 
constrained by the collective properties of the satellite popu- 
lation. 

3.2. The missing satellite problem 

The number of visible structures in our Galaxy constrains 
the stellar mass that CDM subhaloes can form before being 
accreted onto the host and becoming "satellites". It is useful 
to define the star formation efficiency as the fraction of the 
virial mass that a DM halo converts into stars normalized to 
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the universal baryon fraction, i.e. 

Because collision-less CDM simulations do not contain 
baryons, F^, is typically "measured" by associating the most 
massive subhaloes to the most luminous satellite galaxies. 
Due to the divergent nature of the subhalo mass function, 
/j(Mvii-) oc M~iJ.-^, a strong suppression of star formation in 
low-mass haloes is needed in order to reconcile the (small) 
number of visible satellites with the overwhelming number of 
substructures predicted by CDM (a.k.a the "missing satellite 
problem"; Moore et al. 1999; Klypin et al. 1999). E.g. 

in haloes with Mvh < lO'^M© (Kravtsov 2010; see also 
Tollerud et al. 2008 and Koposov et al. 2009). Recent studies 
show that eq. (8) is also consistent with the abundance of lu- 
minous structures found in external galaxies (e.g. Behroozi et 
al. 2010). 

From eq. (6) the formation of a DM core with > O.Irs 
requires 

^■....<-«)a"'-'(:^)"(2|)"' 

Comparison with eq. (8) shows that /v.coi-e/^V.sat ~ 
(Mvi,-/lO"'M0)'/3 < 1 at Mvir < IO'^Mq, implying that the 
formation of DM cores in dSphs leads to an over-abundance 
of luminous satellites (see also Fig. 2). 

However, this comparison has an important caveat, viz. 
that halo masses and concentrations are calculated at z = 0, 
whilst abundance-matching techniques return star formation 
efficiencies at the time of accretion (zacc)- For MW-like haloes 
the accretion redshift of the surviving satellite population 
peaks at Za^c ~ 1 (Zentner & Bullock 2003; Peiiarrubia & Ben- 
son 2005). 

Let us first inspect the energetics of the core/cusp transfor- 
mation as a function of redshift. Because the Universe, and 
therefore the structures that form within, was denser in the 
past, AE tends to increase as we go to higher redshifts at 
a fixed halo mass. This dependence is approximately linear 
during the matter-dominated era (z > 1), because the virial 
over-densities Syn- = pvu/ Pm- 1 — 187r^ were essentially inde- 
pendent of redshift, and because the concentration evolves as 
c(z) « c(z = 0)(1 +z)-' (Wechsler et al. 2002). Hence at a fixed 
halo mass we find that /v,core(z) ~ ^V.corefe = 0)(1 +z) at z > 1. 

Unfortunately, the redshift at which the profile transforma- 
tion takes place (denoted here as Zcoie) is unknown. Given that 
/v.sat is measured at the time of accretion (Zacc), there remains 
the possibility that DM cusps were removed at Zcore ^ Zacc, so 
that Mvir(Zcoi-e) ^ ^vir(-Zacc)- Bccause the energy that baryons 
must supply scales as /S.E oc M^j^., at sufficiently high redshifts 
the formation of cored DM profiles becomes less energetically 
demanding, even though DM haloes were on average denser. 
Below we inspect this possibiUty in detail. 




2 4 6 8 10 



Fig. 3. — Minimum stellar mass required to form a DM core with Tc = \ 
kpc, M*_core, against the stellar mass tuned to reproduce the number of visible 
substructures in CDM simulations, M*,sat. Solid lines show both quantities 
calculated at an accretion time Zacc = 1, which roughly corresponds to the time 
at which the majority of the surviving satellites are accreted onto MW-like 
haloes. Dotted and short-dashed Unes show how the stellar masses compare 
if the cusp removal is shifted to earlier times in the mass evolution of the 
satellites (zcoic = 6 and 10, respectively). Thin and thick lines adopt eoM = 
0.05 and 0.40, respectively. 

3.3. Transformation redshift 

Fig. 3 compares the minimum stellar mass required to form 
a DM core with = 1 kpc, M^^ core, against the stellar mass 
tuned to reproduce the number of visible substructures in 
CDM simulations, My,,sat (eq. 8). 

Models that reconcile the 'missing satellite' and 'core/cusp' 
problems by suppression star formation obey the condition 
^*,core ^ Af*.sat and fall below the horizontal dotted line in 
Fig. 3. Notice that if DM cusps are removed by the time satel- 
lites are accreted (Zcore = Zacc ~ 1 ; solid line) the presence of 
DM cores should be limited to galaxies with > IO^Mq 
for EDM = 0.40 (thick lines), in agreement with the hydro- 
dynamical simulations of G12. 

The threshold between cored and cuspy DM profiles moves 
to fainter luminosities if the transformation is shifted to ear- 
lier times (see §3.2). How early in the past can it be shifted? 
According to Pontzen & Governato (2012) and Teyssier et al. 
(2012) cusps are removed on a time-scale of ~ 1 Gyr, which 
roughly corresponds to the time between multiple in- and out- 
flows of gas during subsequent SNell explosions. Thus in 
a ACDM cosmogony Zcore ^ 6 appears the earliest plausible 
time for the transformation, whereas z ^ 10 corresponds to 
the earliest epoch of star formation in low-mass haloes (Bovy 
& Dvorkin 2012). 

Dotted and short-dashed lines in Fig. 3 show the effects of 
shifting the cusp removal at Zcore = 6 and 10, respectively. In 
order to derive conservative estimates we assume no further 
star formation after the DM core is formed. Halo masses are 
evolved back in time using Wechsler et al. (2002) formalism. 
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These authors show that (i) CDM haloes grow exponentially, 
so thatMvir(acore) ~ A/vir(flacc)exp[-2ac(aacc/flcore- 1)], where 
a = (1 is the scale factor; and (ii) the "collapse-time" of 
the halo (oc) correlates tightly with the virial concentration, 
so that fle ^ 4.1aacc/c(Zacc) independent of the redshift when 
the halo is observed. 

Fig. 3 shows that although shifting the transformation to 
Zcore ~ 6 helps to accommodate DM cores in the bright dSphs 
(i.e. M^, > IQ^Mq), the tension cannot be completely elim- 
inated but is simply shifted to lower luminosities. Notice 
also that the formation of cored halo profiles in dSphs (M^, < 
10^ Mq) requires eoM ~ esn 1- 

4. DISCUSSION 

Our estimates highlight a tension between CDM predic- 
tions and observations on small galactic scales. On the one 
hand, we find that the cored density profile measured in two 
of the brightest MW dSphs points toward an efficient conver- 
sion of primordial gas into stars on the mass scales of dSphs 
(M-jii- < lO'^Mo). On the other, a strong suppression of star 
formation is required on the same mass scales in order to ac- 
commodate the small number of visible satellites with the 
halo mass function predicted by collisionless CDM simula- 
tions. 

It is useful to identify conditions whereby this tension may 
be eased and discuss their compatibility with observations, (i) 
The removal of DM cusps requires a very efficient transfer 
of supernova energy into DM particles (i.e. edm ~ esn)- The 
results of Pontzen & Governato (2012) suggest that the re- 
quired efficiency may be reached if starbursts are sufficiently 
frequent and energetic, (ii) The supernova energy coupling 
(esn) must be order unity. However, such a strong coupling 
may be incompatible with the observed luminosity-metallicty 
relationship of dSphs (Revaz & Jablonka 2012) (iii) Cusps 
must be removed at high redshifts. The cored profiles found 
in Sculptor and Fornax suggest Zcok ^ 6. However, stellar 



ages are at odds with this scenario, as both galaxies show ex- 
tended periods of star formation continuing for 6-7 Gyr in 
Sculptor (de Boer et al. 2012a), and for 12-13 Gyr in Fornax 
(del Pino et al. 2011; de Boer et al. 2012b). Also, even if 
DM cores formed when their haloes had tiny masses, subse- 
quent hierarchical accretion of unperturbed subhaloes would 
tend to re-grow DM cusps (Dehnen 2005). (iv) A possibility 
not explored here is the case of a top-heavy IMF. A higher 
fraction of massive stars would increase the feedback energy 
reservoir without aggravating the over-abundance of luminous 
substructures. 

An alternative route out of the tension is for the satellite lu- 
minosity function to be substantially reshaped by mechanisms 
other than suppression of galaxy formation. Recent hydro- 
dynamical simulations show that baryonic feedback lowers 
the number of surviving satellites with respect to the DM-only 
case (Zolotov et al. 2012). This difference is partially due to 
the enhanced mass-loss rate of satellites embedded in cored 
DM haloes (Penarrubia et al. 2010). Baryonic feedback may 
therefore ease the constraints discussed above by increasing 
the tidal disruption rate of massive (bright) subhaloes moving 
on eccentric orbits. The fact that the star formation efficien- 
cies of the bright MW satellites lie above those dictated by 
abundance matching arguments (Sawala et al. 2012; Boylan- 
Kolchin et al. 2012) may point in this direction. 

Our results also provide compelling reasons to push mea- 
surements of halo mass profiles to galaxies with < 1O*M0, 
as the tension between 'core/cusp' and 'missing satellites' 
problems only increases toward lower luminosity. 
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